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We compute the three point correlation functions for primordial scalar and tensor 
fluctuations in single field inflationary models. We obtain explicit expressions in the slow 
roll limit where the answer is given terms of the two usual slow roll parameters. In a 
particular limit the three point functions are determined completely by the tilt of the 
spectrum of the two point functions. We also make some remarks on the relation of this 
computation to dS/CFT and AdS/CFT. We emphasize that (A)dS/CFT can be viewed 
as a statement about the wavefunction of the universe. 
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1. Introduction and summary of results 

Single field inflationary models predict to a good approximation a Gaussian spectrum 
of primordial fluctuations. The size of non-gaussian corrections is expected to be small 
and was estimated in [§ § . 

In this paper we will compute the correction to the Gaussian answer to leading order 
in the slow roll parameters but with the precise numerical coefficient as well as momentum 
dependence. In single field inflationary models one considers the action of a single scalar 
field coupled to gravity. This action is expanded around a spatially homogeneous solution. 
The leading order term in the expansion is quadratic in the small fluctuations around the 
homogeneous answer. Since a free field is a collection of harmonic oscillators and these 
harmonic oscillators start their life in the ground state one finds that the fluctuations are 
gaussian to leading order. The non-Gaussian effects come from the cubic interaction terms 
in the full action. These interaction terms arise from the non-linearities of the Einstein 
action as well as from non-linearities in the potential for a scalar field. We compute the 
cubic terms in the lagrangian. These cubic terms lead to a change both in the ground state 
of the quantum field as well as non-linearities in the evolution. These two effects can be 
computed in a simple way by following the usual rules of quantum field theory, assuming 
the standard choice of vacuum for an interacting field. 

We parameterize the scalar fluctuations in terms of £ which is the gauge invariant 
variable that remains constant outside the horizon ||. We schematically denote by 7 the 
tensor (or gravity wave) fluctuations. In the slow roll approximation we obtain 



4 (i-i) 

ivi p i 

where e is a slow roll parameter and M.i are homogeneous functions of the momenta of 
degree k~ 6 whose explicit form we give below. The dependence on H/M p i is due to the 
fact that we are looking at the cubic term in the action. Of course, the power of k~ 6 in 
Aii comes from approximate scale invariance. 
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In the limit that one of the momenta in (|1 . 1| ) is much smaller than the other two there 
is a simple argument which determines the three point functions. This simple argument 
can also be used to understand the factors of slow roll parameters in ( |1 . 1| ) . The argument 
is the following. Consider the limit k± <C ^2,3- Suppose the small momentum corresponds 
to a scalar fluctuation £. The fluctuation is frozen by the time the other two momenta 
cross the horizon. So its only effect is to rescale the other two momenta so that we get 
a contribution proportional to the violation in scale invariance of the two point function 
of the two fluctuations with large momenta. In other words the first and third line of 
( p. . 1|) are proportional to the tilt of the scalar and tensor fluctuations respectively (times 
the product of the corresponding two point functions). More explicitly, for the three £ 
correlator we 

where n s is the tilt of the scalar spectrum defined by ~ fc- 3 + n s so that n s is the 
deviation from scale invariance. 

In order to understand the behavior when ki are all of the same order of magnitude 
we need to do the computation by expanding the action. Then the answer is a more 
complicated function of ki but the size of the correlation function does not numerically 
change much. The other two correlation functions in (|1 . 1|) can also be understood in the 
limit that one of the ki is very small through a similar simple argument which we give in 
section 4. 

Another way of presenting the argument is as follows. Since the wavefunction of 
gravity in a space that is approximately de- Sitter is supposed to have the properties of 
a conformal field theory Q0, the three point functions that we computed above can be 
related to correlation functions of the stress tensor in the hypothetical dual CFT. In the 
limit that one of the ki is much smaller than the other two the form of the three point 
function is determined in terms of the two point function by the following argument. If 
one of the ki is very small we can approximate it by zero, so that the corresponding 
insertion of the trace of the stress tensor represents the effects of an infinitesimal rescaling 
of coordinates. So this three point function is determined by how the two point function 

1 We are dropping some factors of (2tt) 3 5(^2 k). These are more explicitly written later. 
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behaves as we rescale the coordinates. This is why the three point function (C 3 ) is equal 
to the tilt of the spectrum of the two point function in the regime fci <C £2,3- 

Komatsu and Spergel have performed an analysis of the detectability of non-gaussian 
features of the temperature fluctuations ||[f§. Their analysis was made for an expected 
signal which had a slightly different k dependence from the one in M.\ above. This probably 
would not change their answer too much. Ignoring this point, one would conclude from 
their analysis that this level of non-gaussianity is not detectable from CMB measurements 
alone. A more explicit discussion is given below. In some models with more than one field 



non-gaussianity can be large JTIJ . 

Finally we point out that these computations can also be used in investigations of 
AdS/CFT and dS/CFT. These dualities can be viewed as a statement about the wave- 
function of the universe. We relate explicitly the computation of stress tensor correlators 
in the dS and AdS case. They are related by a simple analytic continuation. We also 
clarify the relation between stress tensor correlators and the spectrum of fluctuations of 
metric perturbations. 

This paper is organized as follows. In section two we review the standard results that 
follow from the quadratic approximation and give the gaussian answer. In section three 
we expand the action to third order. In section four we compute the three point functions. 
In section five we make some remarks on the relationship of these computations to the 
dS/CFT and AdS/CFT correspondences. 



2. Review of the quadratic computation 

The computation of primordial fluctuations that arise in inflationary models was first 
discussed in |llj] Jl2| |L3| [|14| || jl5[] and was nicely reviewed in fl6[ . 



The starting point is the Lagrangian of gravity and a scalar field which has the general 



form 



12.1) 



S=\J ^g[R-(V<f>) 2 -2Vm 

up to field redefinitions. We have set M~ t 2 = 8ttGn = the dependence on Gat is easily 
reintroduced. 



2 Note that this definition of M p i is different from the definition that some other authors use 
(including Planck). 
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The homogeneous solution has the form 



ds 2 = -dt 2 + e 2p{t) dxidxi = e 2p (-dr] 2 + dxidxi) (2.2) 

where 77 is conformal time. The scalar field is a function of time only, p and <fi obey the 
equations 

3p 2 J-ft + V(<t>) 

p = -\? ( 2 -3) 

=4> + 3p<P + V'(<j)) 

The Hubble parameter is H = p. The third equation follows from the first two. We will 
make frequent use of these equations. 

If the slow roll parameters are small we will have a period of accelerated expansion. 
The slow roll parameters are defined as 



1 fMpiV'V 1 </> 2 1 



2 V V J 2p 2 Mi 

, .. . P (2.4) 

where the approximate relations hold when the slow roll parameters are small. 

We now consider small fluctuations around the solution (|2.3|) . We expect to have three 
physical propagating degrees of freedom, two from gravity and one from the scalar field. 
The scalar field mixes with other components of the metric which are also scalars under 
SO (2) (the little group that leaves k fixed). There are four scalar modes of the metric 
which are <5<yoo, Sga, dgoi ~ d{B and Sgij ~ didjH where B and H are arbitrary functions. 
Together with a small fluctuation, 5(f), in the scalar field these total five scalar modes. The 
action ( |2.1|) has gauge invariances coming from reparametrization invariance. These can be 
linearized for small fluctuations. The scalar modes are acted upon by two gauge invariances, 
time reparametrizations and spatial reparametrizations of the form x l — > x % + e l (£, x) with 
e % = die. Other coordinate transformations act on the vector modest. Gauge invariance 
removes two of the five functions. The constraints in the action remove two others so that 
we are left with one degree of freedom. 
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There are no propagating vector modes for this Lagrangian (2^1). They are removed by gauge 



invariance and the constraints. Vector modes are present when more fields are included. 
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In order to proceed it is convenient to work in the ADM formalism. We write the 
metric as 

ds 2 = -N 2 dt 2 + h lJ (dx i + N l dt)(dx J + N J dt) (2.5) 



2.6) 



and the action (|2.1|) becomes 
S = i J NR( 3) - 2NV + N~ 1 (E i jE ij - E 2 ) + N' 1 ^ - N^^) 2 - Nh ij d i( j)d 

Where 

E ij = l(h ij -V i N j -V j N i ) 

z (2.7) 

E = E\ 

Note that the extrinsic curvature is Kij = N~ 1 Eij. In the computations we do below it is 
often convenient to separate the traceless and the trace part of E^j. 

In the ADM formulation spatial coordinate reparametrizations are an explicit symme- 
try while time reparametrizations are not so obviously a symmetry. The ADM formalism 
is designed so that one can think of h{j and <p as the dynamical variables and N and N l as 
Lagrange multipliers. We will choose a gauge for hij and <p that will fix time and spatial 
reparametrizations. A convenient gauge is 

5<f> = , = e 2p [(l + 2C)<% + 7«] , 9 <7y = , lu = (2.8) 

where ( and 7 are first order quantities. £ and 7 are the physical degrees of freedom. £ 
parameterizes the scalar fluctuations and 7 the tensor fluctuations. The gauge ( |2.8| ) fixes 
the gauge completely at nonzero momentum. In order to find the action for these degrees of 
freedom we just solve for N and N l through their equations of motion and plug the result 
back in the action. This procedure gives the correct answer since iV and N l are Lagrange 
multipliers. The gauge ( |2.8j ) is very similar to Coulomb gauge in electrodynamics where 
we set diAi = 0, solve for Aq through its equation of motion and plug this back in the 
actionS. 

The equation of motion for N l and iV are the the momentum and hamiltonian con- 
straints 

V l [N-\E)-5)E)} = Q 

(2.9) 

#(3) _ 2V - N- 2 {E i:j E lJ - E 2 ) - N~ 2 (j) 2 = 



4 As in electrodynamics in Coulomb gauge we will often find expressions which are not local 
in the spatial directions. In the linearized theory it is possible to define local gauge invariant 
observables where these non-local terms disappear. 
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where we have used that 8<p = from ( j2.8|) . We can solve these equations to first order by 
setting N l = ditp + where diN l T = and N = 1 + N x . We find 

AT X = 4 , JV*=0, V = - e - 2 ^ + x , d 2 X = ^C (2-10) 

P P 2p 2 



In order to find the quadratic action for £ we can replace ( |2.10[ ) in the action and 
expand the action to second order. For this purpose it is not necessary to compute N or 
N l to second order. The reason is that the second order term in N will be multiplying 
the hamiltonian constraint, ^ evaluated to zeroth order which vanishes since the zeroth 
order solution obeys the equations of motion. There is a similar argument for N l . Direct 
replacement in the action gives, up to second order, 

S =i I e p+c (l + i)[-48 2 C - 2(d() 2 - 2Ve 2p+2C ] + 

I . . (2-11) 

+ 

where we have neglected a total derivative which is linear in ijj. After integrating by parts 
some of the terms and using the background equations of motion ( |2.3| ) we find the final 
expression to second order! 

S= \ I dtd 3 x^[e 3p ( 2 - e p (d() 2 } (2.12) 
^ J p 

No slow roll approximation was made in deriving (|2.11|) . Note that naively the action 



( |2.11| ) contains terms of the order £ 2 , while the final expression contains only terms of the 
form e£ 2 , so that the action is suppressed by a slow roll parameter. The reason is that 
the C fluctuation would be a pure gauge mode in de-Sitter space and it gets a non-trivial 
action only to the extent that the slow roll parameter is non-zero. So the leading order 
terms in slow roll in (|2.11|) cancel leaving only the terms in ( |2.12j ). A simple argument for 
the dependence of ( |2.12| ) on the slow roll parameters is given below. 

Since ( |2.12|) is describing a free field we just have a collection of harmonic oscillators. 
More precisely we expand 



In order to compare this to the expression in [16| set v = —zQ in (10.73) of [16|. 
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Each Ck(t) is a harmonic oscillator with time dependent mass and spring constants. 
The quantization is straightforward jl7| . We pick two independent classical solutions £^ (A) 
and (%. l *(t) of the equations of motion of ( |2.12| ) 

Then we write 

( l ;(t)=C k l (t)at + (t(t)a_ i; (2.15) 

where a and are some operators. Demanding that and a obey the standard creation 
and annihilation commutation relations we get a normalization condition for £j£ . Different 
choices of solutions are different choices of vacua for the scalar field. The comoving wave- 
length of each mode A c ~ 1/k stays constant but the physical wavelength changes in time. 
For early times the ratio of the physical wavelength to the Hubble scale is very small and 
the mode feels it is in almost flat space. We can then use the WKB approximation to solve 
( p. 14] ) and choose the usual vacuum in Minkowski space. When the physical wavelength 
is much longer than the Hubble scale 

X phys H = ^- > 1 (2.16) 

the solutions of ( |2.14| ) go rapidly to a constant. 

A useful example to keep in mind is that of a massless scalar field / in de-Sitter space. 
In that case the action is S = | J H~ 2 r]~ 2 [{d Tt f) 2 — (df) 2 ] and the normalized classical 
solution, analogous to , corresponding to the standard Bunch Davies vacuum is JHJ 

fS = ^|=(1 - ikv)e ikr > (2.17) 

where we are using conformal time which runs from (— oo,0). Very late times correspond 
to small \r]\ and we clearly see from (|2.17|) that f cl goes to a constant. Any solution, 
including ( p. 171) , approaches a constant at late times as rf ~ e _2p , which is exponentially 
fast is physical time. In de-Sitter space we can easily compute the two point function for 
this scalar field and obtain! 



Wv)f*{v)) =(^) 3 S 3 (k + k')\f c k l (v)\ 2 = (2n) 3 5 3 (k + k')^(l + k 2 V 2 ] 

tt2 

~(2n) 3 5 3 (k + k')—r for krj < 1 
2k 6 



(2.18) 



6 In coordinate space the result for late times is (f(x,t)f(x',t))sim — -^—z\og(\x — x'\/L) 
where is an IR cutoff which is unimportant when we compute differences in / as we do in actual 
experiments. 
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We now go back to the inflationary computation. If one knew the classical solution to 
the equation ( |2.14| ) the result for the correlation function of £ can be simply computed as 

<Cje(f)C*(t)> = (2n) 3 5 3 (k + k')\Q l (t)f (2.19) 

If the slow roll parameters are small when the comoving scale k crosses the horizon 
then it is possible to estimate the late time behavior of (|2.19| ) by the corresponding result 



in de-Sitter space ( |2.18| ) with a Hubble constant that is the Hubble constant at the moment 
of horizon crossing. The reason is that at late times £ is constant while at early times the 
field is in the vacuum and its wavefunction is accurately given by the WKB approximation. 



Since the action (|2.12|) also contains a factor of (p/p we also have to set its value to the 
value at horizon crossing, this factor only appears in normalizing the classical solution. In 
other words, near horizon crossing we set f = where / is a canonically normalized field 
in de-Sitter space. This produces the well known result 

<wc*(«)> ~ ( 2 -) 3 ^+^i|§ ( 2 - 2 °) 

where the star means that it is evaluated at the time of horizon crossing, i.e. at time 
such that 

p{U)e p{u) ~ k . (2.21) 



The dependence of ( 12.20]) on t* leads to additional momentum dependence. It is conven- 



tional to parameterize this dependence by saying that the total correlation function has 
the form k~ 3+ris where 

n " = *4 log (|) ~ ~JT log (|) =-2(-^- + h = 2(77 " 3e) (2.22) 

As it has been often discussed, after horizon crossing the mode becomes classical, in 
the sense that the commutator [£, Q — > exponentially fast. So for measurements which 
only involve £ or £ we can treat the mode as a classical variable. 

After the end of inflation the field (p ceases to determine the dynamics of the universe 
and we eventually go over to the usual hot big bang phase. It is possible to prove M |T6[] that 



£ remains constant outside the horizon as long as no entropy perturbations are generated 
and a certain condition on the off-diagonal components of the spatial stress tensor is 
obeyedB. These conditions are obeyed if the universe is described by a single fluid or by a 



7 The condition is didj(8Tij — hSijTu) = . 
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single scalar field. We should mention that for a general fluid the variable £ can be defined 
in terms of the three metric as above ( |2.8| ) in the comoving gauge where Tf = oi. In the 
case of a scalar field this implies that 8(f) = 0. This gauge is convenient conceptually since 
the variable £ is directly a function appearing in the metric. We see that the variable £ tells 
us how much the spatial directions have expanded in the comoving gauge, so that to linear 
order ( determines the curvature of the spatial slices = Ak 2 X ||19|| . This variable £ is 



very useful in order to continue through the end of inflation since it is defined throughout 
the evolution and it is constant outside the horizon. An intuitive way to understand why £ 
is constant is to note that the conditions stated above imply that two observers separated 
by some distance see the universe undergoing precisely the same history. Outside the 
horizon (where we can set k = in all equations) C is just a rescaling of coordinates and 
this rescaling is a symmetry of the equations. 

Other gauges can be more convenient in order to do computations in the slow roll 
approximation. A gauge that is particularly convenient is 

8(f> = <p(t, x) , hij = e 2p (8 ij + jij) , dciij = , 7« = (2.23) 

where we have denoted the small fluctuation of the scalar field by (p. In order to avoid 
confusion, from now on (f> will denote the background value of the scalar field and <p will be 
its deviation from the background value. We expect that in this gauge the action will be 
approximately the action of a massless scalar field ip to leading order in slow roll. Indeed, 
we can check that the first order expressions for iV and N l are 

= K = S iX , * X =^±U V ) ,2.24) 



where the <p subindex reminds us that Aq^, AT* are computed in the gauge (|2.23| ). We see 



that these expressions are subleading in slow roll compared to (p. So in order to compute 
the quadratic action to lowest order in slow roll it is enough to consider just the (V</?) 2 
term in the action (|2.1| ) since V" is also of higher order in slow roll. This is just the 



For readers who are familiar with Bardeen's classic paper |lq| , we should mention that the 
gauge invariant definition of £ is £ = h + (J^T x b! — A)TL 2 /{Ti 2 — 7€) where Ti = p' and primes 
indicate derivatives with respect to conformal time and h = Hl + Ht/3 with A, Hl, Ht defined 



in [18 1 . In circumstances where £ is conserved then it also reduces to the definition in terms of 
Bardeen potentials in ||, |16| | (actually Chere = — Cthere)- The gauge choice that makes h = £ is 
= or, using the equations of motion, h = pA. 
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action of a massless scalar field in the zeroth order background. We can compute the 
fluctuations in ip in the slow roll approximation and we find a result similar to that of a 
scalar field in de-Sitter space ( |2.18|) where the Hubble scale is evaluated at horizon crossing. 
After horizon crossing we can evaluate the gauge invariant quantity £. This is most easily 
done by changing the gauge to the gauge where (p = 0. This can be achieved by a time 
reparametrization of the form t = t + T with 

T = -% (2.25) 

where t is the time in the gauge ( |2.8|) and t is the time in (|2.23| ). After the gauge trans- 
formation ( |2.25|) , we find that the metric in ( |2.23|) becomes of the form in ( |2.8| ) with 



C = pT = -£ip (2.26) 



Incidentally, this implies that x i n ( P-24j ) is the same as x i n (|2-10| )- So the correlation 



function for ( can be computed as the correlation function for ip times the factor in ( |2.26|) . 
In order to get a result as accurate as possible we should perform the gauge transformation 
( |2.26|) , just after crossing the horizon so that the factor in (|2.26| ), is evaluated at horizon 



crossing leading finally to (|2.20|) . In principle we could compute ( from ip at any time. If 
we were to choose to do it a long time after horizon crossing we would need to take into 
account that <p changes outside the horizon. This would require evaluating the action ([2.1] ) 



to higher order in the slow roll parameters. Of course, the dependence for ip outside the 



horizon is such that it precisely cancels the time dependence of the factor in (|2.26|) so that 
( is constant. 

In summary, the computation is technically simplest if we start with the gauge ( |2.23| ) 
and we compute the two point function of tp after horizon exit and at that time compute the 
C variable which then remains constant. On the other hand the computation in the gauge 
Q2.8|) is conceptually simpler since the whole computation always involves the variable of 



interest which is (. In other words, the gauge (|2.23| ) is more useful before and during 



horizon crossing while the gauge (2.8) is more useful after horizon crossing. 



These last few paragraphs are basically simple argument presented in |T3|]. The com- 
putation of fluctuations of <p in de-Sitter produces fluctuations of the order ip = and 
then this leads to a delay in the evolution by St = —<p/p (see ( |2.25|) ) which in turn gives an 
additional expansion of the universe by a factor £ = p5t = — £(p. This additional expansion 
is evaluated at horizon crossing in order to minimize the error in the approximation. 



11 



We now summarize the discussion of gravitational waves [pDJ . Inserting (2^8) in the 
action and focusing on terms quadratic in 7 gives 

S = I j ^ PAiijAiii ~ ePdl ^ d nn\ ( 2 -27) 

As usual we can expand 7 in plane waves with definite polarization tensors 

^=/(0E4(%I(*)^ (2-28) 

where en = k l €ij = and efJk)efJk) = 25 ss r. So we see that for each polarization 
mode we have essentially the equation of motion of a massless scalar field. As in our 
previous discussion, the solutions become constant after crossing the horizon. Computing 
the correlator just after horizon crossing we get 

(7|7,<> = (2nfS\k + k')± ||<W (2.29) 

pi 

where we reinstated the M v \ dependence. We can similarly define the tilt of the gravita- 
tional wave spectrum by saying that the correlation function scales as k~ 3+nt where rtt is 
given by 

n t = k-^\ogpl = -% = -2e (2.30) 

UK 



3. Cubic terms in the Lagrangian 

In this section we compute the cubic terms in the lagrangian in two different gauges. 
We do this as a check of our computations. The first gauge is similar to ( |2.8| ), which is 
conceptually simpler since one works from the very beginning with the ( variable in terms 
of which one wants to compute the answer. We need to fix the gauge to second order in 
small fluctuations. We achieve this by setting to zero the fluctuations in <p and we writing 
the 3-metric as 

5<j) = 

1 (3-1) 

hij = e 2p+2C hij , det h = 1, = (<% + 7^ + — TiiTij H ) 

where 7^ = d^ij = to second order. The term proportional to 7 2 was introduced with 
the purpose of simplifying some formulasi. Note that it is necessary to define hij only to 
second order since any third order term in hij will not contribute to the action. 

9 We can define the gauge condition as ft (log h)ij = 0. 
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The second gauge that we choose is 



(f) = (f)(t) + tp(t,x) 

1 (3-2) 
hij = e 2p hij , det h = 1 , = (<% + + -^lulij H ) 

again with 7^ = <9i7ij = 0. In appendix A we work out in detail the change of gauge which 
gives the relation between the C, lij variables and the p : 'jij variables. Here we summarize 



that discussion. We denote by t the time in the gauge (|3.2|) and by t the time in the 
gauge (|3.1] ). t and t are related by a time reparametrization of the form t = t + T(t,x). 
The function T should be such that from Sep 7^ in ( |3.2| ) we get Sep = in ( |3.1| ). This 
determines 

Starting from ( |3.2| ) this time reparametrization produces a new metric. A spatial 
reparametrization then carries it to the form in (|3.1|) to second order. After all these 
steps we find the following relation between the 99 and ( variables (for 7 = 0) 

C =pT + \pT 2 - \d l Td J Te- 2p + \d iX diT+ 

where T is defined in ( |3.3|) and x is defined in ( |2.24| ). This expression will be useful for 
comparing results computed in different gauges. The change of variables relating 7^ and 
7ij can be found in appendix A. 



3.1. Evolution outside the horizon to all orders 

It is possible to show that ( and 7 (defined in (|3.1|) ) are constant outside the horizon. 
For this purpose it is enough to expand the action to first order in derivatives of the 
fields but to all orders in powers of the fields. We will assume that = 1 + 5 N where 
8N has an expansion in derivatives that starts with a first order term. Similarly we will 
assume that N l is of zeroth order in derivatives, that that VjNi are of first order in 
derivatives. These assumptions are consistent with the structure of the hamiltonian and 
momentum constraints (|2.9| ). We can then expand the Hamiltonian constraint to first 
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order in derivatives and solve for SN to first order in derivatives and all orders in powers 
of the fields^ 

2VSN = 2p(3C - V,iV l ) (3.5) 

We can now evaluate the action. On a solution of the Hamiltonian constraint the action 
reads 



S = J VhN(R 3 -2V)= J Vh(-2V - 2V5N) = 
= J e 3p+ ^(-6p 2 + <j> 2 - 6Cp) = -2 / d t (e 3 ^p) 



(3.6) 



where we neglected the term involving because it is of second order in derivatives, we 
used (|3.5|) , we integrated by parts the term involving ViN 1 and we used (p73|). Therefore 
(fT6|) is a total derivative and can be ignored. 

In conclusion, £ and 7 are constant outside the horizon. The reason is that outside 
the horizon we can neglect all spatial derivatives. Since we also showed that the expansion 
in powers of time derivatives starts at second order we conclude that constant £ and 7 are 
solutions of the equations of motion to all orders in powers of £, 7 outside the horizon. 
This fact was derived in a different way in Of course, the intuitive explanation of 

this fact is clear, after exiting the horizon all regions evolve in the same fashion, the only 
difference between them is how much one has expanded relative to the other. This fact 
makes it clear that the definition of £ in (|3.1| ) |PT | is the correct non-linear generalization 
of the variable introduced in ||. 

3.2. Cubic term for three scalars 

We now expand the action up to cubic order in £. It turns out that it is only necessary 
to know N or N l up to first order. The third order terms in N , N l would be multiplying 
the constraints evaluated at zeroth order. The second order terms in N, N l multiply the 
constraints evaluated to first order, which vanish due to the first order expressions for iV 
and N l . Up to total derivatives in time and space we find 
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S=J e^(l + ^)[-2d 2 C - (d() 2 ] + e 3 ' +3 ^C 2 (l - i)+ 



+ C 3P+3C 



hdid^did^ - (dV) 2 )(l - -) - 2d^d l Cd 2 ^ 
I p 



(3.7) 
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Note that if we expand this to linear order in the fields it only agrees with our linearized 



expressions 
11 



2.10) to first order in derivatives. 



[21| call a our (. 
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where we expand the exponentials so that we keep only terms of up to third order in £, 
and ip is defined in ( |2.10| ) and is of first order in £. 

Something that is not obvious from ( pi.7| ) is that the effective cubic interaction term 
is of second order in slow roll. Schematically it is of the form e 2 ^ 2 ^, while the action 
( P-7| ) appears to be of order e . Factors of e are very important since they will determine 
whether this effect is measurable or not [Q. 

An easy way to see that the effective action is of order e 2 is to compute the cubic term 
in the action in the gauge ( |3.2|) , which leads to 



30 3 3 4> 5 3 4>V" 3 V" 3 ^ 2 . <f 2o2 ,, v 

+ j7(-<pd i d j xd i d j x + (pd\d 2 x) 



Only the terms in the first line of ( |3.8| ) contribute to leading order in the slow roll ap- 
proximation. The term proportional to V" was considered in but we see that it is 
subleading in the slow roll approximation. By using the first order relation between £ and 
(p ( |2.26|) we see that the first line of (|3.8| ) leads to an effective action which is of order e 2 
in the £ variables. On the other hand, in the action ( |3.8|) , it is not obvious that there is 
any variable which stays constant outside the horizon. Indeed there are <p s couplings that 
typically lead to evolution of the perturbations outside the horizon 0. 

Since one property (the constancy of Q is clear in one gauge while the other (the 
fact that the interaction is of order e 2 ) is more clear in the other the reader might have 
some doubts about one or both statements. In order to dissipate all doubts about these 
statements we start from the cubic term in £ in ( |3.7|) , do a lot of integrations by parts and 
drop total derivative terms to obtain 



+ /(C)| 



1(f) 1(f) 2 



+ l^e 3 Pd i d jX d i d jX ( (3.9) 
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where the first line indicates the leading order terms, which are of order e 2 as expected. 
Note that x is °f order e, see ( |2.10p . The second line indicates higher order terms in the 



slow roll expansion. Finally the third line denotes terms that are proportional to the first 
order equations of motion ( |2.14j) . They can be removed by performing a field redefinition 
of the form 

C =Cn - /(Cn) 



1 6 .o 1 



/,2 



2 6p 4 p l 



p 4 p z 4 p z 2 p 2 p 

(3.10) 

where we have written the explicit expression for /0. After this field redefinition the 
action in terms of ( n is just the first two lines of ( |3.9|) . 

By comparing the field redefinition ( ]3.1(J| ) with (|3.3|), (|3~^|), ( |2-26j ), we find that 



Cn = ~^ (3.11) 

with no quadratic correction. This provides a consistency check on our computations since 
it is clear that the actions for ( n and ip have the same form to leading order in slow roll 
(the first line of (|3]8|) should be compared with the first line of ( |3.9|) ) and indeed the ( and 
( n are related as we expect by the corresponding change of gauge. The agreement between 
the two forms of the action should persist to all others in slow roll but we did not verify it 
explicitly. 



It should be noted that the field redefinition (|3.10|) does indeed matter for our com- 
putation since we are interested in computing expectation values of ( and not of Q n . The 
reason is that ( is the variable that stays constant outside the horizon while ( n does not. 
This last fact follows from the fact that £ is constant and the form of ( |3.10| ) where some 
of the coefficients of the quadratic terms are time dependent. We can also see from the 
action ( |3.9| ) that the second line involves a term with only one time derivative on ( n . This 
term leads to evolution of ( n outside the horizon consistent with what is expected from 
( |3.1C)|) . Note that only the terms in the middle line of (|3.10|) are important outside the 
horizon. 



12 



Note that it does not matter whether set C or Q n in the quadratic terms. 
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An easy way to perform the computation is to compute using the <p or ( n variables up 
to a few Hubble times outside horizon exit of the relevant modes and then change variables 
to the £ variables where it is clear that there is no evolution outside the horizon. 

In order to perform the computation of the three point function we will use a variable 
Cc defined through 



C = C + -^C 2 + -^C 2 + -^V 2 (C d 2 C) + ■ 



(3.12) 



where the dots indicate terms that vanish outside the horizon or are higher order in the 
slow roll parameters. In terms of £c the action becomes 



Ss = J P 5p K 2 c d- 2 C c + 



(3.13) 



where the dots again indicate terms of higher order in slow roll. The last term in ( |3.12j ) 
comes from terms proportional to the equations of motion that arise when we integrate by 
parts the first line in (|3.9f) in order to get (|3.13|) . 



3. 3. Cubic term for two scalars and a graviton 
We start with the computation in the gauge 



S 



- e^ijdiCdjC - le 3p (3( - -)",,<),<),,.- + le^dmjdid^drt (3.14) 



Again, it is easiest to understand the dependence of the action on the slow roll pa- 
rameter by computing the action in the gauge ( |3.2| ) , where the leading contribution comes 
from 

S = J )pijdi(pdj(p H (3.15) 

where the dots indicate terms that are of higher order in slow roll. We can then conclude 
that, despite appearances, ( |3.14[ ) should be of order e. 

Indeed, after some integrations by parts, we find that (|3.14|) becomes 



S= J l^ lij d i (d i (+ 



+ ^P d ^. jdixdjX + ^ e ZPj. jd . Cd . x + /( C> 7 ) 



SL 



SL 



(3.16) 
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where the first line indicates the leading order term in slow roll, which indeed has the same 
slow roll dependence as ( |3.15| ), once (|3.11| ) is taken into account. The second line contains 



a higher order term in the slow roll approximation as well as well as terms proportional to 
the equations of motion. These terms can be removed by field redefinitions which we give 
explicitly in appendix A. These have the form that we expect from changing the gauge from 
(|3.1|) to (|3.2|) . It turns out that these field redefinitions are not important after horizon 
crossing and hence are not important for our computation. 

3.4- Cubic term for two gravitons and a scalar 

Let us first do the computation in the gauge (|3.1| ) • By direct substitution in the action 
and after some integrations by parts and dropping total derivatives we get 

S = [ ^[e^^^ + e^da^da^-^^da^x 

J P £T (3.17) 

-Cnj-c — + ••• 

As usual the second line can be removed by a field redefinition. This field redefinition is the 
same one that we have to use to go from the gauge ( |3.ip to the gauge ( |3.2| ) , as is discussed 
in more detail in appendix A. When we do the computation in (|3.2|) we get directly the 
first line of ( |3.17|) , after taking into account ( 3.11 ). It is curious that the form of the first 



line in ( |3.17| ) is rather similar to the first line in (|3.9| ). As we did in that case, in order to 
perform the computation it is convenient to do the further field redefinition 

C = Cc " + T6 a_2 ^ a2 ^) + ' ' • ( 3 - 18 ) 

where the dots indicate terms that vanish outside the horizon. Then the action becomes 

S = J \^pe 5p i^jd- 2 (c + ■■■ (3.19) 
where the dots indicates terms that are higher order in the slow roll approximation. 

3. 5. Cubic term for three gravitons 

The computation of the term involving three gravitons is the same in the gauge ( |3.1| ) 
or the gauge ( |3.2| ), since after setting the scalar to zero we are changing the metric in 
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precisely the same way. We note that the only terms in the action that contribute come 
from 

S=\J e 2 ?{R + E)El) (3.20) 

This has the same form as the result we would have obtained if we had done the compu- 
tation in flat space, except for the factor of e 2p . The variable hij was defined in terms of 
7 in such a way that there is no cubic term involving two time derivatives. This implies 
that when we integrate by parts we will not need to use the equations of motion and 
therefore there will not be any field redefinitions. In flat space we know from the form of 
the scattering amplitudes that we can reduce the effective vertex to a term involving only 
spatial derivatives. We give its explicit form in the next section. 



4. Computation of three point functions 

In this section we compute the three point functions using each of the interaction 
lagrangians that we found above. Before describing the computation let us make a couple 
of general remarks. 

First let us notice that we are computing an expectation value and not a transition 
amplitude. We want to compute, in the interaction picture, 

<C 3 (*)> = (Ur n lcHt)U int (t,to)) , U int = T e -^>" t(t ' )dt ' (4.1) 

where to is some early timeEl, and T denotes a time ordering. We have suppressed the 
spatial dependence. To first order this gives 

<C 3 W> = -* f 'dt'([C 3 (t),H int (t')}) (4-2) 

Jto 

For the cubic terms Hi nt = —Li nt after we remove all terms proportional to the equations 
of motion by a field redefinition. 

Our second point is to note that we want to compute the expectation value in the 
vacuum of the interacting theory, not the vacuum of the free theory. When we do computa- 
tions in Minkowski space we also have to take this into account. This can be automatically 

13 So we are not computing (T£ 3 e %nt ) which is what we compute when we have scat- 

tering amplitudes in mind. When we compute scattering amplitudes field redefinitions do not 
change the answer. In our computation they do. 
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taken into account by deforming the t' integration contour so that it includes some evo- 
lution in euclidean time which projects on to the true vacuum. Fortunately in this case 
we can apply a similar procedure to select the vacuum. The basic reason is that at early 
times the physical wavelength is very small and we feel in Minkowski space, so we want 
the vacuum for these high energy modes to be what it is in the interacting theory in this 
approximately Minkowski space. In de-Sitter space this is the Hartle Hawking prescription 
for the vacuum . In practice this will translate into a choice of contour for the integral 



in (|4.2|) .lli The evaluation of the integral in Q4.2| ) reduces to an integration of the cubic 
action evaluated on the classical solutions of ( |2.14|) Ei. Since we do not have the solutions 
for a general potential it is useful to choose a method that minimizes the errors in the ap- 
proximate evaluation of the integrals. These errors are minimized if we split the integrals 
in ( |4.2j ) as an integral over the region outside the horizon, the region around horizon cross- 
ing and the region deep inside the horizon. In the last region the fields oscillate rapidly 
and after our continuation to Euclidean space there is no contribution. In the region near 
horizon crossing we approximate the solutions by those of de-Sitter space ( |2.17| ) and we 
use the action in the form that shows the leading slow roll dependence, such as in ( |3.13| ), 
( |3.16|) or ( |3.19| ) . After we exit the horizon we know that £ and 7 are constant so we switch 
to those variables. This is taken into account by the field redefinitions we talked about. 
Then in the region well outside the horizon the fields are constant and the integral (|4.2| ) 
vanishes when we do the computation in the £, 7 variables. Below we proceed with this 
computation for the various cases. Due to momentum conservation there are basically two 
distinct kinematic configurations, the k{ can be all of the same order of magnitude of one 
of the ki is much smaller than the other two. We will consider these two cases separately. 

4-1. Three scalar s correlator 

Note that if we have a field redefinition of the schematic form C = Cc + M,c then the 
correlation function will contain two terms 

(C(x 1 )C(x 2 )C(x 3 )) = (Cc(zi)Cc(x 2 )C c (x 3 )> + 2A [(ax 1 )C(x 2 ))(ax 1 )C(x 3 )) + cyclic] (4.3) 



14 Other choices of vacua in de-Sitter space were discussed in [p3| . In inflation the admixture 



of these other vacua is expected to be small [24| and the leading contribution to the three point 
function comes from the usual vacuum. 

15 This is true only after performing field redefinitions to eliminate terms proportional to the 
equations of motion. Otherwise we need to take into account total derivative terms that arise 



when we go from (3/7) to (|3.9| ), for example. 
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The first term is computed by ( |3.13| ) . The second comes from the field redefinition ( |3.12|) . 
By performing different field redefinitions we can reshuffle the contributions between these 
two terms. Let us first compute the contribution from the action ( |3.13| ). As we explained 
above we evaluate this term in de-Sitter space with parameters corresponding to those of 
horizon exit. In de-Sitter space the contribution of an action of the form ( |3.13| ) has the 
form 

(CcCcCc)=(27r) 3 5 3 (V^)TT7^T3T4 i / dr)k\k\j k ^ + permutations + c.c. = 

LL( 2k i)(pi J-oo 

=(27rf8 3 h) |U kj , with k t = k l+ k 2 + h 

(4.4) 



Note that fcj = \ki\. The choice of integration contour in (|4.4|) is such that the oscillating 
piece in the exponent becomes exponentially decreasing. In other words we change r\ — > 
r\ + ie\r}\ for large \rj\. This choice of contour is the one corresponding to the standard 
vacuum of the interacting theory. 

After adding the contribution of the field redefinitions we get the final result for the 
three point function 

(C^Cfe) = ygg^A (4.5) 

where the star indicates evaluation at horizon crossing and 



4>*p* i P* 



(4.6) 



In writing ( |4.5| ) and (^.6| ) we have assumed that all /c's are of the same order of magnitude 
so that the moment of horizon crossing does not differ too much between the different 
modes. Due to momentum conservation the other possibility is that one of the /c's is much 
smaller than the other two and these last two would be of the same order of magnitude. 
So we consider the configuration k% k\ ~ k 2 - The mode labeled by k^ crosses the 
horizon much earlier than the other modes. By the time that fci 2 cross the horizon £3 is 
constant. The only effect of the £3 fluctuation will be to make the comoving scales k\p 
cross the horizon at a slightly earlier time 5t* = —Qz/p*. This will produce a change in 
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the fluctuations with momenta due to the time dependence of the slow roll factors in 
( |2.20|) . In conclusion we obtain 

^-n s 4Ck 3 C_^)(Ck t Ck 2 ) (4.7) 

where now * indicates the moment that ki, h% cross the horizon and *' indicates the time 
when ^3 crosses the horizon (which is earlier). In the second line we point out explicitly 
that this two point function involves the tilt of the scalar spectrum n s * (|2.22| ), evaluated 
at t*. The prime in the first two lines of ( |4.7|) means that we omit the factor (27r) 3 5(0). It 
can be checked that (|4.5| ) , with (|4.6|) goes over to ( |4.7|) in the overlapping region of validity 
which is when ks is small but not so small to change the slow roll parameters appreciably. 
The first two lines in (|4.7| ) are valid to all orders in slow roll parameters in the regime 

Our result Q4.5| ) is of the same order of magnitude as in [|J , but the k dependence 
as well as the precise numerical coefficients are different. The reason is that Q considers 
only effects due to non-linear evolution but do not consider the change in vacuum. Both of 
these effects are of the same order of magnitude so they should be both included and are 
intimately linked. Our result obeys consistency condition explained in the above paragraph 
while that in ||] does not, presumably because not all relevant effects were included. 

Spergel and Komatsu [§0[j2l| did an extensive analysis of the measurability of the 
three point function. They assumed that the three point function had the form that would 
follow from a field redefinition of the formll! 

C = C 9 - \fNLC 2 g (4-8) 

where ( g is gaussian. Their analysis can be roughly summarized by saying that this would 
be measurable if Jnl is bigger than around 5. This constraint comes mainly from cosmic 
variance if we assume that we measure the CMB up to the angular scales that the Planck 
satellite will measure them. See PS[ for a detailed discussion of this point. Our final result 



16 Spergel and Komatsu defined = <E> g + /nl^- The factor of 5/3 arises in the relation 
between the gauge invariant newtonian potential $ and £ during matter domination, £ = — 1$. 
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does not have the momentum dependence that would follow from (|4.8|) . In order for 
that to be the case we would need that all terms in ( f4.6| ) were proportional to £V kf which 
is clearly not the case. So we cannot recast our computation as a computation of /jvl- 
Nevertheless we can define a /c-dependent /jvl as 



A 



NL 



3(4E 4 *f) 



5 

12 



+ ^(2 + f(k)) 

pi 



.-(n s + f(k)n t ) (4.9) 



where f(k) has a range of values < / < f . f(k) is a function of the shape of the triangle 
made by ki and it goes to zero when two sides become much larger than the third and it 
becomes 5/6 when the ki form an equilateral triangle. 

Assuming that this k dependence of the three point function does not significantly 
change the analysis in IH1 1|25|1 , we unfortunately conclude that it will not be possible 



to see this effect purely from the CMB. Actually, the discussion of [25] makes sense for 
for > 1 where we can neglect the non-linearities in the gravitational evolution after 
horizon reentry, some of these effects were discussed by [p6l|[271||2"8|] . In other words, to 
measure Jnl < 1, one has to include the leading non-linear effects in the whole evolution 
until we measure the temperature of the CMB. 



4-2. Two scalars and a graviton correlator 

This computation is rather similar to the one we did above so we will not repeat all 
the details. Let us note that there is no field redefinition that is important at late times 
so that we only need to evaluate the integral that arises from the interaction term in the 
first line of ( [Up . 

This gives 



i Pi p 2 



* ^elklkUl 



(4.10) 



where the transverse and traceless polarization tensor is normalized to efj-e^- 
/ is 

I =Re 



f - ik lT j)(l - ik 2V )(l - iks^e^ 

J — oo V 



I = -k t + 



Y.i>jkikj kik 2 k 



25 ss ' and 



(4.11) 



+ 
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The integral in ( 4.11 ) diverges at r\ — > but the divergence is purely imaginary so that / 
is finite with our choice of contourtH. 

The dependence on the slow roll parameters is such that the three £ correlator and 
the 7C 2 correlator are of the same order of magnitude. After horizon reentry the amplitude 
of the gravitational waves decays so that for high I we still expect the three £ correlator 
to dominate. 

Let us now consider the correlation function in the limit k\ <C &2, &3- When &3 
are about to cross the horizon the gravity wave with momentum k\ is already frozen so 
that the fluctuations of £ will be those that we expect in this deformed geometry. The 
main effect of the deformation is to change k 2 —>■ k 2 — ^ijk 1 ^ in the correlation function 
of two (s. This reasoning leads to 

<^ 2 c* 3 >' - - (^7: Sl >^^(c 5j cs s > 



\ ) \/_^ "01.5 01.3,2^2 M 2 . 9 2 2 



where * denotes the time when /c2,3 cross the horizon while *' denotes the time when k\ 
crosses the horizon. 

We see that ( |4.12| ) and ( 4.10 ) are consistent in the overlapping region of validity. This 



is a consistency check of the computation. 
4-3. Two gravitons and a scalar correlator 



The evaluation of this correlator using (|3.18|) and Q3.19|) is very similar to the one of 
the three scalar correlator. We obtain 

/a s 2 S3 \ _ (2?r) 5(J2 k i) f> r 1,3 1 , / ,2 , ;2n , M l J 2 « 3 (A 1^ 

<vwv _ n(2^) 3 m p v 4 /ei + 2 /ei(/C2+/e3j+4 Je ^ (4 - ldj 

In the case that k\ <C /c2,3 we also find that the correlation function is given by the 
derivative of the slow roll factor in the correlation function of two tensor fluctuations. We 
get 

~ -MC^C^/'frgTg} (4 . 14) 
~(fcr) ^^-^^-^^ 

17 In order to evaluate this integral it is convenient to note that Re[— i J° oo dr]r]~ 2 (l — iw)e lv ] = 
— 1 with our contour prescription. 
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where we have used that k^ ~ k% so that e||e*| ~ 25 S2S3 . The * indicates horizon crossing 
for &2, &3 and the *' indicates horizon crossing for k\ which happens earlier. In the second 
line we emphasized the dependence of this three point function on the tilt of the gravity 
wave spectrum. We see that ( |4.14|) agrees with fl4.13|) in the overlapping region of validity. 



4-4- Three graviton correlator 

The three graviton correlator is a very similar computation. The algebra involving 
polarization tensors is the same as in flat space so that we can use the flat space result. 
We will need to do the same integral as in ( |4.11|) . The final result is 

<-0£t£> = ( 2 -)^(E^^ n ^(-4)(^^^^^/0^ (4-15) 

where / is given in ( |4.11|) , and tijk is given by the flat space formula (see for example J29[) 

Uji = kiSji + k J 3 5u + k[d l3 (4.16) 

We can compute this in the limit k\ -C /c2,3 in a way similar to what we did for the 
case of a graviton and two scalars 

(7^7^) = (27r) 3 £ 3 (y ^ ^f' i.^i.-d'fc'H (4.17) 
which indeed agrees with ( [4 . 1 5| ) in the overlapping region of validity. 



5. Remarks on AdS/CFT and dS/CFT 

5.1. AdS/CFT 

The computation that we did above was done with inflation in mind, but the same 
mathematical structure arises if one considers a single scalar field with a negative potential. 
In the slow roll case, the background will be a slightly deformed anti-de-Sitter space. This 
can be understood as a slightly deformed conformal field theory. In other words, a non- 
conformal field theory which is almost conformal. An incomplete list of references where 
situations of this sort were considered is f30|j ||31|| [p2[ ||33|| []34]] ||35|| . Here we just mention a 



few results that are relevant for us, for a review see [^6|. The variables 7 s that we used 
above are associated to the traceless components of the stress tensor while the variable 
( is associated to the trace of the stress tensor. More precisely, we have a coupling of 
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the form / -^[2C_^(k) + 2 7 * -T s (£)], where T s is defined by an expression similar 
to ( |2.28[ ), with 7 — > T. The fact that the definition of the scalar mode depends on the 
gauge is translated into the fact that in a field theory with a scale we can either change 
the dimensionfull coupling constant or we can change the overall scale in the metric. It is 
common to fix the coupling and change the metric, which then relates Q to the trace of 
the stress tensor. Alternatively we can fix the metric and change the coupling constant. 
In the field theory we do not have two independent operators, we have only one operator 
related by the equation 

2Ti = (3 x O (5.1) 

where f3\ is the beta function for the coupling A which appears in the field theory La- 
grangian in front of the non-marginal operator as J XO. The operator O is the one coupling 
to <fi and the operator 2T? couples to £. The factor of slow roll that relates the correlators 
of C and 4> is precisely the factor (3\ appearing above [ p7| . 

From the computations in the previous sections we can also compute the correlation 
function of stress tensors and trace of the stress tensor in non-conformal theories. Depend- 
ing on whether the slow roll approximation is valid or not we would need to use different 
formulae in those sections. 

Two point functions of the trace of the stress tensor were considered in the AdS 
context in [|33| [^2| |34] [ 35| . The derivation of the effective action for the corresponding field 
in AdS identical to the one in the dS context. Similarly, computations of three point 
functions in AdS can be done by performing minor modifications to the above formulae. 
We will be more explicit below. 

Now we will review the AdS/± computation (see ||38|| for a review) so that we can 
contrast it clearly to the dS± computation. 

Let us consider a canonically normalized scalar field in Euclidean anti- Sitter space 
(EAdS^) which is the same as hyperbolic space. The action is 

S = R\ dS J d ~[{d z f) 2 + {dfY] (5.2) 

In order to do computations it will be necessary to consider classical solutions which go to 
zero for large z and obey prescribed boundary conditions at z = z c . In momentum space 
these are 

/s = % (1 + kz c )e-^ ' k = ]k] (5 ' 3) 
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where /S is the boundary condition we impose at z = z c . One should then compute the 



action for this solution as a function of the boundary conditions. Inserting (|5.3| ) into (|5.2j 
integrating by parts and using the equations of motion we get 



_ , d?k^ 1 2 o 1 df n 



(2tt) 3 2 AdS J -k J k Zc (l + kz c 



f 1 2 f0 r^_,3 

J (2^2 HAdsf -^ [ z c k + 



(5.4) 



where the dots indicate terms of higher order in z c . The term divergent in z c is local in 
position spaced and it is viewed as a divergence in the CFT which should be subtracted 
by a local counterterm. The term independent of z c is non-local and gives rise to the two 
point function 



(o(k)o(k')) EAdS i) " z 



(2nfS(k + k')R 2 AdS k 3 (5.5) 



/°=o 



Where Z is the partition function of the Euclidean CFT which is approximated by Z 
e~ Sc \ with S in fl57§). 

5.2. dS-CFT 



The dS/CFT was proposed 0[| in analogy with AdS/CFT H@g|F|. The dS/CFT 



postulates that the wavefunction of a universe which is asymptotically de-Sitter space can 
be computed in terms of a conformal field theory. More precisely, we have the formula 

n 9 } = Z[g] (5.6) 

where the left hand side is the wavefunction of the universe for given three metric and the 
right hand side is the partition function of some dual conformal field theory. Actually the 
left hand side has rapidly oscillating pieces which can be expressed as local functions of the 
metric. We discard these pieces since they have the interpretation of local counterterms in 
the CFT. Here we are thinking of de-Sitter in flat slices (or Poincare coordinates) and we 
are imagining that all fields start in their life in the Bunch-Davies vacuum. This determines 
the wavefunction at least in the context of perturbation theory. If we were considering 
global de-Sitter space then our discussion would be valid in a small patch in the future 



18 It is proportional to £ / dx 3 ±(dff. 
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where it can be approximated by the Poincare patch and the memory of the particular 
state that could have come from the far past is lostM This point of view follows simply 
from the discussion in || in analogy with the standard discussion in Euclidean AdS where 
the same formula 05.61) is valid0. Nobody has found a concrete example of this duality and 
there are some suspicions that such a duality should not exist [45| . All we will do here is 
to do some computations on the gravity side in order to get some insight on the properties 
that this hypothetical CFT should have. If an example were found, then it would be a 
more powerful way of computing the wavefunction that semiclassical physics in de-Sitter 
or nearly de-Sitter space. Note that an observer living in eternal de-Sitter space will not 
be able to measure two point correlators such as ( |2.29| ) or the wavefunction ( |5.6| ) which 
involves distances much larger than the Hubble scale. Only so called "metaobservers" can 
measure these . On the other hand if the universe is approximately de-Sitter for a while 
and then inflation ends and we go over to a radiation or matter dominated universe then 
these correlation functions become observable. In fact, we are metaobservers of the early 
inflationary epoch [46|| . 

In [f?l ||47|1 the relation between CFT operators and fields in the bulk was explored 
and various ways of defining operators were considered. It was found that given a scalar 
field in the bulk one could define two operators with two conformal dimensions differing 
by A_|_ — A_ = d where d is the dimension of the CFT. If the field we are considering 
in the bulk is the metric then it is clear that the corresponding operator is the stress 
tensor and it should have dimension d. Indeed we will see that this agrees precisely with 
what we expect from the prescription ( |5.6|) . Below we explain more precisely how this 
computation is related both to the inflationary computation ( |2.2y| ) and the corresponding 
EAdS computation. 

The first step is to compute the wavefunction as a function of a small fluctuation in a 
massless scalar field /. Since / is a free field, which is a collection of harmonic oscillators, 
all we need to do is to compute the wavefunction for these harmonic oscillators. We want 



19 The information of the state coming from the asymptotic past in global dS is contained on 
modes whose angular momenta, /, on the sphere is fixed, assuming the evolution is non-singular 
and in the context of perturbation theory. On the other hand, we are focusing on modes with 
I 3> 1 when we look at the Poincare patch. 

20 In AdS / CFT formula ( |5.6| ) arises in the Euclidean context when we think of Euclidean time 
as the direction perpendicular to the boundary. ^ can then be interpreted also as the Hartle- 



Hawking wave function [22]. See [42] [43] [44 for more on this point of view. 
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to compute the Schroedinger picture wavefunction at some time r\ c as a function of the 
amplitude of the field /. The wavefunction is given by a sum over all paths ending with 
amplitude / and starting at the appropriate vacuum state. Since the action is quadratic 
this sum reduces to evaluating the action on the appropriate classical solution. We choose 
the standard Euclidean (Bunch-Davies) vacuum for the fields at early times. The classical 
solution obeying the appropriate boundary conditions is 

o (1 - ikrj)j kr > 

k (1 _ ikr] c )e ikr i- 1 ' ' 

The boundary conditions at large r\ are the ones that correspond to the statement that 
the oscillator is in its ground state, which can be defined adiabatically at early times. 
The condition is that the field should behaves as e lkr> for \r)\ — ► oo. Note that /_£ 7^ /i 
since the boundary condition we are imposing at early times is not a real condition on 
the field f(rj,x) 0. This is one of the many ways to think about the harmonic oscillator 
wavefunction. When we evaluate the classical action on this solution we get 

d*k 1 k 2 ( } 



(2tt) 3 2 a * 1 rj c 1J ~ kJk 

Note that we are dropping an oscillatory piece at I77I —>■ 00 which is equivalent to slightly 
changing the contour of integration by r\ — > rj + ie. This is the standard prescription for 
the vacuum state of a harmonic oscillator. 
Notice that under 

r] = iz , R d s = iRAdS (5.9) 

the formulas Q5j7p and Q go into (^) and (pH ). The fact that O goes into (^) is 
intimately related to the statement that when the mode has short wavelength it is in the 
adiabatic vacuum. A consequence of this fact is that the two point function computed using 
dS^ differs by a sign from the corresponding one in Euclidean Adsj^. More explicitly we 
have 

d 2 Z 



(0(k)0(k f )) 



dS^ 



~(27r) 6 8(k + k')R 2 ds (-k*) (5.10) 
/°=o 



21 There is nothing wrong in considering a complex solution since all we are doing is to evaluate 
a functional integral by a saddle point approximation. 

22 In other dimensions there are extra is that appears in the relation. 
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We can easily check that this is the analytically continued version of (5J3) under 

Now let us understand the relation between the wavefunction computed in ( |5.8| ) , which 
is ^ ~ e %Scl and the expectation values that appeared in our earlier discussion (|2.18j ). Of 
course, the relation is that (f 2 ) = f Vf f 2 \^(f)\ 2 . We see that only the real piece in iS 
contributes. This has a finite limit at late times. The divergent pieces in (|5.8|) are all 
imaginary and do not contribute to the expectation value. The functional integration over 
/ gives again (|2.18| ). There is a crucial factor of 2 that comes from the square of the 
wavefunction, so that the relation between (|2.18|) and (|5.10| ) is not a Legendre transform. 

Our previous discussion focused on a scalar field and its corresponding operator O. All 
that we have said above translates very simply for the traceless part of the metric and the 
traceless part of the stress tensor, since at the linearized level the action for the graviton 
in the traceless transverse gauge reduces to the action of a scalar field ( |2.27| )( |2~T28| ). We 
are defining the stress tensor operator as 

Tij[x) s «*M = *»M (5.11) 

which is a standard definition for a Euclidean field theory.il In this case the divergent 
term in ( |5.8| ) can be rewritten as — i-^- J d 3 x\/hR^ 3 \ Note that there is a factor of i. We 
want to remove this by a counterterm in the action of the Euclidean CFT. These factors 
of i are related to the fact that the renormalization group transformation in the CFT 
should be appropriately unitary since this RG transformation corresponds, in the context 
of perturbation theory, to unitary evolution of the wavefunction in the bulk. If we define 
the central charge of the CFT in terms of the two point function of the stress tensor we get 
a negative answer. This negative answer has a simple qualitative explanation. We know 
that the wavefunction in terms of small fluctuations is bounded, in the sense that it is of 
the form e _a '^' with a positive, since each mode is a harmonic oscillator with positive 
frequency. This sign implies a negative sign for the two point function of the stress tensor. 
Similarly the trace of the stress tensor is related to the derivative of the wavefunction with 
respect to (. 



23 One might want to define it with an i so that Tji = i Sj*gj~jl ■ This definition might be natural 
given that the counterterms (which represent the leading dependence of the wavefunction) are 
purely imaginary. In any case, it is trivial to go between both definitions. 



30 



After we understood the relation between two point functions of operators and expec- 
tation values of the corresponding fluctuations we can similarly understand the relation 
between three point functions. The wavefunction has the form 



^ =Exp 



]- I d 3 xd 3 x'(0(x)0(x'))f(x)f(x') + 

, , i (5-12) 

- / d 3 xd 3 x'd 3 x"(0(x)0(x')0(x"))f(x)f(x')f(x") 



6 

where we emphasized that derivatives of \& give correlation functions for the corresponding 
operators. The expectation values in momentum space are related by 

(ftf-k) 



2Re(O u O_ n )' 

2Re(Or Or Or V ^ A3 > 

i f f f y — ' fc i fc 2 k 3 i 

\JkJk 2 Jk 3 / -Yl.(-2Re(0 % p_ k )>) 

where the prime means that we dropped a factor of (27r) 3 5(^ k). And Re indicates the 
real part. The factors of two come from the fact that we are squaring the wavefunction 
( |5.12|) . Notice that this explains why (TT) ~ c while (77) ~ 1/c where c ~ —R^gM^. 

Now consider three point functions. For example, consider the three point function of 
the traceless part of the stress tensor. This can be computed directly in dS^ by inserting 
the classical solutions (|5.7[) into the cubic terms in the action. This gives 



M 2 1 

< T £ T £ T £> = ^f^h)^ (-^)i^^t ijlti , j/v )i (5.i4) 



where / is defined in (|4.11[ ) . The result in EAdS^ is the same as above except for a minus 
sign, which can be understood as coming from (|5.9|) . When we perform this computation 
we need to drop a local divergent term which is proportional to ^VhR^ . We did not 
have any divergence in ( f4.17| ) due to the fact that we were computing the square of the 
wavefunction while in ( |5.14j ) we are computing the third derivative of the wavefunction. 



Of course, we can compute directly (|5.14|) from ( [4.17| ) using (|5.13|) . So in order to compute 



three point functions of the stress tensor in the hypothetical three dimensional field theory 
corresponding to a nearly dS± spacetime all we need to do is apply formula ( |5.13|) to our 
results in section four. To go to the corresponding expectation values in EAdS^ we just 
need to multiply all dS^ results by a minus sign which comes from R^g — > —R\dS an< ^ 
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all correlators of the stress tensor have such a factor in front in the tree level gravity 
approximation. 

Some of the points we explained above are specific to the four dimensional (IS4 case. 
The situation in dS§ is rather interesting. The computation of fluctuations for a massless 
scalar field gives, outside the horizon, 

(f % f %l ) ~ H 3 (27r)*6(k + %)~ , H = (5.15) 

On the other hand the wavefunction ^ ~ e lS has the form 

iS = --R ds [ -^fy\{^-h 4 log(- Vc k)-i-k 4 + ak 4 } (5.16) 
2 ds J (27r) 4 k ~ k l 2r% 4 K ' 8 J 

where a is a real number. Note that the only term contributing to ( |5.15|) is the real term 
proportional to k 4 . All other terms are purely imaginary. From ( |5.16| ) we can compute 
the non-local contribution to the two point function which gives 

(0(k)0(k')) dS5 ~ (2TT) 4 5(k + k')iRl±k 4 \ogk (5.17) 

The EAdS$ answer is given by the analytic continuation ( |5.9| ). Notice that the i is due 
to the fact that we have an odd number of powers of R d s and is consistent with the fact 
that the logarithmic term in the wavefunction is purely imaginary. For the stress tensor 
this gives an imaginary central charge and imaginary three point functions. It is rather 
interesting that the two point function ( |5.15| ) is related to a local term in the wavefunction, 
namely the term proportional to k 4 , which is the only real term. In other words, the non- 
local piece in the wavefunction which determines the stress tensor seems unrelated to the 
local piece which determines the expectation value of the fluctuations. In other words, 
dS^/CFT^ would tell us how to compute the non-local piece in the wavefunction but will 
give us no information on the local piece. On the other hand from the inflationary point 
of view we would be interested in computing ( |5.15| ) which depends on the local part of 
the wavefunction, or the partition function of the CFT. Maybe in dS/CFT we are only 
allowed to use imaginary counterterms, then the field theory should be such that it allows 
the computation of the finite real local parts in the effective action. Note that the real term 
in ( |5.16|) arises in the analytic continuation ( |5.9|) from the term in the EAdS$ wavefunction 
that is proportional to k 4 \og(z c k) — > — \k 4 + k 4 log(— t] c k). So still, in some sense, the 
real part of the wavefunction (|5.16| ) is intimately related to the non-local term in the 
wavefunction. It looks like this will be the situation in all odd dimensional dS spaces. 
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The AdSz case studied in M seems special because there is no bulk propagating graviton. 



Stress tensor correlators in dS/CFT were also studied in [f48flp7 



Now let us reexamine the three point functions of stress tensor operators in the limit 
that one of the momenta is much smaller than the other two. We can then approximate 
the small momentum by zero. This zero momentum insertion of the stress tensor can be 
viewed as coming from an infinitesimal coordinate transformation. So we then know that 
the three point function is going to be given by the change of the two point function by 
this coordinate transformation. For example, an insertion of the trace of the stress tensor 
at zero momentum is equivalent to performing a rescaling of the coordinates without 
rescaling the mass scale of the theory. Then the three point function will be given by the 
scale dependence of the two point function. In other words 

(2Tf(0)O(k)O(k')) = -k*^-(0(k)0(k')) (5.18) 

This is the reason why three point functions in this limit are proportional to the tilt of 
the scalar and tensor spectra respectively, see ([O]) (|4.14|) . There is a similar argument for 
the insertion of the traceless part of the stress tensor at zero momentum. Formula (|5.18| )is 
valid to all orders in slow roll. 

Notice that in order to compute observable quantities from dS/CFT we will need to 
square the wavefunction and integrate over some range of values of the couplings and the 
metric of the space where the CFT is defined. In other words, in order to compute some 
physically interesting quantity it is not enough to consider the CFT on a fixed 3-manifold 
but over a range of three manifolds. This is the reason that expectation values in dS 
are not simply given by analytic continuation of the ones in EAdS ffl] even though the 
wavefunction and correlation functions of operators are given by analytic continuation.!! 
This makes it clear that even if dS/CFT is true there is no causality problem, one is 
not fixing the final state of the universe. One fixes it as an auxiliary step in order to 
compute the wavefunction but in order to compute probabilities we need to sum over all 
final boundary conditions. A slightly different integral over boundary conditions arises 
also in the EAdS context when we consider certain relevant operators [^9| , or double trace 
operators ||50||. In those cases this integration is the same as a change in the boundary 



24 This analytic continuation is very clear for fields with 2mR < is d < d. For fields with mass 
above this bound it is not so clear what the right prescription is. In this paper we focus our 
attention on the easy case. 
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condition. Note that this is not what happens in the dS context since we have the square 
of the wavefunction. One might conjecture that dS expectation values are given by two 
CFTs (one for \& and one for coupled together in some fashion. Note that then it is 
not clear if we should view the resulting object as a local field theory since in the resulting 
object is not defined on a fixed manifold since in order to compute expectation values we 
need to integrate over the three metric. The two copies of the CFT that we are talking 
about arise just at the future boundary, so these two copies are different than the two 
copies talked about in [j6| ||47|1 1|48|| . In global coordinates in addition we have the past 
boundary. Throughout this paper we have ignored the past boundary since we focused 
on distances larger than the Hubble scale but smaller than the total size of the spatial 
slice. In the Hartle and Hawking prescription for the wavefunction of the universe the past 
and future parts of the wavefunctions are complex conjugates of each other since the total 
wavefunction is real ]2^]. It is natural to suspect that these two pieces can be thought of 
as \1/ and \&* in our discussion above. 
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Appendix A. Second order change of variables between the two gauges 

In this appendix we work out explicitly the change of variables to second order between 
the gauge (|3.1|) and the gauge ( jO]) . Let us denote by i the time coordinate in ( |3.2|) , by t 
the one in ( |3.1| ). The time reparametrization is t = t + T. First let us find the first order 
change of variables between the two gauges. It is easy to see that we do not need to do any 
spatial reparametrization to first order and all we need to do is a time reparametrization. 
The value of T\ at first order, as well as the relation between variables is 

rp Vl Cl J- P /A -f\ 

P <j) 

where the subindex reminds us that it is a first order relation. 
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Now we work this out to second order. In order to go from the gauge ( |3.2|) where ip 
is not zero to ( |3.1] ) where <p is zero we need to do a time reparametrization determined by 
the equation <f>{t + T(t)) + ip(t + T(t)) = (f)(t) which gives, to second order, 

T = _%_\^ + W (A.2) 



Under this change of variables we find that the metric in (|3.2j ) becomes 

ft r = e 2 P (t+T) . + ~ . . {t) + 1 ~ u ~ y + ~ . + ^ , T + ^ a . T _ e -2p d . TdjT) (A 3) 

where we have set N = 1 in some second order terms and is given in ( |2.24| ). This 
metric does not yet obey the gauge condition in (pTl|) . The violation of the gauge 
condition is due to some second order terms, since we already saw that at first order we 
do not need to do a spatial reparametrization. The terms responsible for this violation are 
the last four terms in ( |A.3| ). In order to make it obey those gauge condition it is necessary 
to do a spatial reparametrization where x l = x l +e l (x, t), where e l is of second order. The 
condition that e l should obey is that 

Shlj + die 3 + dj€ l = 2a5ij + % , difiij = , fJ-u = (A.4) 



Where bh\^ represents the last four terms in ( ]A.3| ) . In order to solve this equation for e l it 



is convenient to separate e l = die + e\ where die\ = 0. Taking the trace and didj of (|A.4| ) 
we obtain 

4a =6h r u - d^didjSh^ 

= - d l Td l Te~ 2p + 2Ar;a,T + d- 2 d i d j (d i Td j T)e- 2p (A.5) 
- 2d- 2 d l d J (%d 3 T) - d^didjT 

We can similarly compute what e 1 are and then find from (|A.4|) . It turns out that for 
our purposes we can write 

Hij = 5hij + rest (A. 6) 



where the last terms are given by all other terms in ( A.4 ). These terms vanish when 



integrated against a function which is traceless and divergenceless, which means that these 
extra terms do not contribute to our computation. In conclusion, after making this spatial 
reparametrization we can put the metric in a form such that it obeys the gauge ( |3.1| ) to 
second order. 
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We find then that the final field redefinition is given by 



C =p(t + T) - p{t) + a 

. (A.7) 

where T, a and //^ are given above. 

It is convenient to define a variable ( n through the relation ( 3.11|) . Here we are 



thinking of £ n as a convenient parameterization of the variable p. Then the explicit change 
of variables between the two gauges is 



1 <t> 1 <b 2 „ 2 



2 0p 4 

1 • 1 e~ 2p 1 e _2p 1 1 

+ -CnCn - T^2~(^Cn) 2 + T —^-d~ 2 didj (di( n dj( n ) + --^Xn^Cr, 

p 4 p 2 4 p 2 J J 2 p 



---9 ^didjidiXndjCn) ~ -7-lijdidjCn ^ A ' 8 ^ 
I p 4 p 

1 • e~ 2p 1 

+ ~7«Cn T^OiCndjCn + - (diXndjCn + ^X^Cn) 

f) J f) z f) 

we see that only the first line of each field redefinition is non- vanishing outside the horizon. 
So fortunately we do not need to take all these terms into account in the computations in 
the paper. We did check however that the terms proportional to the equations of motion 
that arise when we integrate by parts the lagrangian in the gauge ( |3.1| ) in order to make it 



look more like the lagrangian in ( |3.2|) are indeed precisely the ones that lead to the above 
field redefinitions. 
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